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Abstract. We solve analytically the ellipsoidally expanding fireball hydrody- 
namics with source terms in the momentum and energy equations, using the 
non-relativistic approximation. We find that energy transport from high p t 
jets of gluons to the medium leads to a transient, exponential inflation of the 
fireballs created in high energy heavy ion collisions. In this transient, inflatory 
period, the slopes of the single particle spectra are exponentially increasing, 
while the HBT radius parameters are exponentially decreasing with time. This 
effect is shown to be similar to the development of the homogeneity of our Uni- 
verse due to an inflatory period. Independently of the initial conditions, and 
the exact value of freeze-out time and temperature, the measurables (single 
particle spectra, the correlation functions, slope parameters, elliptic flow, HBT 
radii and cross terms) become time-independent during the late, non-inflatory 
stages of the expansion, and they satisfy new kind of scaling laws. If the expan- 
sion starts with a transient inflation caused by the gluon wind, it leads naturally 
to large transverse flows as well as to the the simultaneous equality, and scaling 
behaviour of the HBT radius parameters, R S ide ~ Rout ~ Riong ~ tfy/TfJm. 
With certain relativistic corrections, the scaling limit is Tfy^Tf /m t , where rrit 
is the mean transverse mass of the pair. 
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1. Introduction 

Recently, there were a number of QCD based perturbative calculations that indi- 
cated the energy loss of gluon jets with high transverse momentum due to multiple 
scattering during the course of their penetration through hot and dense hadronic 
matter. This phenomena, the jet quenching goes back to papers of Gyulassy, Pliimer 
and Wang [ [j], [|] , recently studied in the papers of Gyulassy, Levai and Vitev (GLV) [ 
pf. It is clear from these works, that the jet quenching mechanism results in the 



2 



T. Csorgo and J. Zimanyi 



energy loss or the depletion of the number of hadrons from the high transverse 
momentum part of the spectrum. Due to the conservation of the total energy, this 
implies that the number of particles with small total momentum has to increase. 
See ref . [ g for the summary of recent developments in this field from the point of 
view of perturbative QCD. 

In this paper, we focus on the phenomenological consequences of energy and 
momentum transfer to the soft, hydrodynamically behaving, small transverse mo- 
mentum part of the single particle spectrum. We present a new family of solutions 
of non-relativistic hydrodynamics with source terms. This family of solutions is 
a generalization of the results of refs. [ ^, ||, |?J to ellipsoidally symmetric, ex- 
panding fireballs that are subject to energy and momentum transfer from the non- 
hydrodynamically behaving modes. A constant rate of energy pumping is shown to 
lead to an exponential inflation of the principal axis of the fireballs, while momen- 
tum transfer leads only to linear increase. These effects are correlated in time to 
the transition of high energy gluons through the nuclear medium, hence, by default, 
they are also transient phenomena, existing only in the initial phase of the expan- 
sion. We point out the similarity of this mechanism to the existence of a transient 
inflatory period in the expansion of the early Universe. 

Inflation of the Universe after the Big Bang leads to a homogeneous and flat 
solution of Einstein's equation. We show here that an inflation of hadronic fireballs, 
caused by the gluon wind in the initial stage of the Little Bangs of heavy ion 
collisions also leads to the simultaneous equality, spherical symmetry and scaling of 
the effective, measurable source sizes (frequently referred to as the HBT radii). This 
effect connects the physics of the smallest observable scales to the physics of the 
largest observable ones, underlining the universal nature of the laws of physics and 
the scale independence of the local conservation laws, governing the hydrodynamics 
of these expansions. 



2. There blows the gluon wind ... 

Consider the non-relativistic hydrodynamical problem, as given by the continuity, 
Euler and energy equations, appended with source terms in the momentum and the 
energy conservation laws: 

d t n + V(vn) = 0, (1) 
<9 t v + (vV)v = -(V(p + pc))/(mn), (2) 
cV + V(ev) +pVv = j G , (3) 

where n denotes the particle number density, v stands for the non-relativistic (NR) 
flow velocity field, e for the energy density, p for the pressure and in the following 
the temperature field is denoted by T. The above set of equations are supplemented 
by the equations of state (EoS). The parameters of the EoS vary in various domains 
of the (T, /i.) plane, corresponding to various phases of matter. A general, and 
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analytically solvable family of equations of state is given by 

p = XpnT - B , e= \ e nT + B. (4) 

The parameters k — A e /A p and B characterize the EoS for a broad variety of ma- 
terials: e.g. a non-relativistic ideal gas yields k = 3/2 and B = 0. Softening of 
the EoS can be modeled by increasing the values of k in certain temperature do- 
mains. A parameter B > can be used to describe a kind of generalized bag model 
equation of state (vacuum pressure) and to describe phase transitions, similarly to 
the class of models considered in refs. [ g|. Hence we present generalizations of the 
results in refs. [ g, ||] and references therein, for a possible description of phase 
transitions during the time evolution of the exploding fireballs. 

The source terms pa and jq in eqs. (§,[§) stand for the contribution of the 
gluon wind that blows from the center of the collision zone through the material 
for a finite period of time. This wind does not carry conserved charges, and we 
assume that the energy and the momentum transfer is proportional to the number 
density n of the fluid, as well as to the average local velocity of the fluid elements. 
We assume that the gluons are proceeding essentially with the speed of light, so the 
relative velocity distribution depends only on the breadth of the local momentum 
distribution, which in turn is characterized with the magnitude of the temperature 
T. We define a model that is analytically solvable by assuming that 

Pg = 3v{t)nT, (5) 
3G = 3E{t)nnT, (6) 

where the coefficients j v and Je depend only on the time variable t. The time 
dependence of the gluon wind can be expressed with the help of the dimensionless 
functions j v (t) and jE(t) in such a way that they are non- vanishing only within 
a finite period of time t\ < t < £2- The equations of the usual non-relativistic 
hydrodynamics are recovered in the case of vanishing source terms, Je(£) = jv{t) = 
0. 

A simple analytic solution of the non relativistic hydrodynamic equations with- 
out source terms (i.e. je = jv = 0) was given for the spherical and cylindrical 
symmetric cases in refs. [ |t| and [ ^| more than 20 years ago, and a recent series 
of papers [ |l2|, [l3| |[ [7]] generalized these solutions to various (Gaussian, donut 
shaped or oscillatory) density and corresponding temperature profiles, for various 
symmetry classes ranging from spherical to ellipsoidally symmetric cases. Ref. 
|l2| also considered the case of Gaussian density profiles, spherically symmetric ex- 
pansions and source terms in the energy equation, and pointed out the existence of 
an inflatory solution with an exponentially increasing radius parameter for the case 
corresponding to js(t) — jo > 0. Here we generalize this solution for the case of 
ellipsoidal symmetry and study in detail the consequences for the observables. 

We choose n, v and T as the independent variables. We find that the NR hydro 
equations with source terms are solved by the following self-similar, ellipsoidally 
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symmetric density, temperature and flow profiles: 
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n(t,r') = n -^exp ( -- % - Jx_ - | , , 7) 
V ' ; l V I 2X 2 2Y 2 2Z 2 i ' w 

v'(*,r') = (jk, | A (8) 



T(t) = T <(*)^J ' ^ 
Ti(t) = /sT [l + j fl (t)]expQf j E (u)du^ (10) 

where the variables are defined in a center of mass frame K' , but with the axes 
pointing to the principial directions of the expansion. This is the frame correspond- 
ing to the System of Ellipsoidal Expansion (SEE) , introduced in ref. [ ^| , where the 
relationship between the observables in SEE and in the CMS of the fireball were 
analyzed in great details. 

The time dependent scale parameters are denoted by (X, Y, Z) the typical vol- 
ume of the expanding system is V — XYZ, and the initial temperature and volume 
are To = T(to) and Vq = V(to), and no is a constant. The time evolution of the 
radius parameters X, Y and Z is equivalent to the classical motion of a particle in 
a non-central, time dependent potential: 

XX = YY = ZZ=^±^y\ (11) 

which correspond to a Hamiltonian motion of a mass point at (X, Y, Z) driven by 
the following non-central, time dependent potential: 

* - Pj± ^ ]1+ Mf)"- <12) 

where {P x , P y , P z ) = m{X, Y, Z). See refs. [ 1 1 for a similar problem, the solution 
of non-relativistic fireball hydrodynamics for self-similarly expanding, ellipsoidally 
symmetric case without source terms in the energy and the Euler equations. 

In all the cases, the normalization constant of the temperature increases in 
time corresponding to reheating during inflation, given by eq. (|Io|). This is due to 
the external sources in the energy and momentum balance equations. This effect 
competes with a decreasing factor, governed by the expansion of the system, the 
rate of which is controlled by the parameter k. Note that for relativistic, massless 
particles, the speed of sound is simply c 2 = 1/k. 

From this Hamiltonian formulation, the effects of the equation of state and the 
effects from the source terms in the energy and momentum balance equations can 
be qualitatively analyzed, and numerical solutions are also very easily obtained. 
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Softening the equation of state decreases slope of the repulsive potential in the cor- 
responding mechanical problem. As the coordinates (X, Y, Z) correspond to the 
characteristic scales of the expanding ellipsoids in the hydrodynamical system, soft- 
ening the equation of state means slowing down the expansion. It is very interesting 
to observe the phenomena of a "stall" during a first order rehadronization transition 
of a quark-gluon plasma [ [2lL |2^, ^3). In our case, the acceleration of the fluid 
does not stop even for the softest equations of state, 1/k — » 0, although the rate of 
expansion slows down substantially, according to eq. (|TT|). 

3. Inflation generates strong transversal flows 

From the analysis of the solution of fireball hydrodynamics with inflation, it follows 
that three different effects may lead to the speeding up of the transverse expansion: 

i) Hardening the equation of state (corresponding to a decrease in the value 
of k) . This does not change the type of the expansion, the coordinates (X, Y, Z) 
increase as power-law type of functions of the time, as the slope of the potential 
term in eq. (|l2|) is increased. However, decreasing k makes the algebraic coefficients 
of the expansion larger. 

ii) Momentum transfer. A positive value of the factor j v (t) is also effective in 
increasing the height of the potential, however, the effect is only linear in j v . 

Hi) Energy transfer. This is the most effective method to inflate the fireball, 
as the increase in the volume changes its form from a characteristic power-law type 
of expansion to an exponential shape. Similar effect was discussed for the case of 
spherically symmetric fireballs in ref . [ [l2| . 

4. Inflation, the gluon wind and scaling of HBT radii 

The observables: the shape, the slope parameters of the single particle spectra, 
the various flow coefficients: first, second and third flows, as well as the various 
HBT radius parameters were calculated analytically in ref. [ |6) for the case when 
the source terms were all set to zero. We found that the introduction of the source 
terms to the energy and the Euler equation does not change the form of the solution, 
only the time evolution of the source parameters is modified. Hence, all the results 
of ref. [ |(| are generalized to the case of inflating expansions. We present some 
numerical examples of the time evolution of the observables, using equations (8-37) 
of ref. [ |(|, that we do not recapitulate here. Although the steps are trivial, the 
results are literally surprising. 

4-1. Inflating hydro solutions with a gluon wind 

In this subsection we present the analysis of the expansion for the initial case of a 
heavy ion reaction, when energy and momentum transfer to the soft modes is sup- 
posedly present, as a consequence of the penetration of the highly energetic, pertur- 
bative gluons through the fireball. For simplicity, we assume that the momentum 
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transfer is negligibly small, j v (t) = 0. We assume that the QCD jet-quenching 
processes and the related energy transport to the soft momentum modes can be 
characterized by an approximately time independent coefficient jE(t) — jo, valid in 
a finite time interval (ti.,^)- In principle, the functions j v (t) and j_e(0 should be 
determined from QCD based Monte-Carlo simulations, if a realistic simulation is at- 
tempted. Here our aim was not yet a detailed fitting of the data, but the qualitative 
understanding of the behaviour of measurablcs, hence we analyzed the behaviour 
of the observables for these reactions. As high pt gluons are emitted from the be- 
ginning of the high energy heavy ion collisions, it is convenient to choose t\ = 0. 
The end of the hard collisions can be estimated by the maximum of the penetration 
time, t 2 ~ 1.12A^' 3 ' /j. For clarity, we investigate analytically the case when the 
fireball is initially spherical, X = Y = Z = R . 

From eq. (^) , we find that the radius of the fireball inflates exponentially, and 
the coefficient of the temperature also increases exponentially, 

R{t) = Roexp^j^-jot^j , (13) 
Ti(t) = T ij0 exp(j i), (14) 

T(t) = m^f^T^^jot). (15) 

Hence the local temperature T(t) increases exponentially for all n > values, be- 
cause the exponentially increasing volume term cannot compensate the even faster, 
exponential increase of the Tj(t) factor. This effect is referred to as reheating during 
the inflation, or reheating caused by the gluon wind. 

Let us consider how such a scenario effects the experimentally measurable slope 
parameters, T e g and the HBT radius parameters, i? c ff- Following the lines of refs. [ 
||, [l^ j, these observables are given as 

T eS = T + mR 2 = (T a + j^-mR 2 )eM^^jot), (16) 

R cS ~ 1 T^exp (--Jl—j t\ , (17) 

jo \3 + 2Ky V tti \ 3 + 2k J 

which implies that the effective slope parameters increase, while the effective ho- 
mogeneity regions decrease exponentially with increasing time. These features of 
the inflatory period remain qualitatively similar for non-spherical initial conditions 
too, as shown in Figs. 1 and 2. 

Hence inflation by the gluon wind may lead to the development of strong trans- 
verse flows as well as small HBT radius parameters. However, after the penetration 
of the nuclei through each other, the production mechanism for high momentum 
gluons comes to an end, and this implies that the inflation by the gluon wind will be 
over. After this time, the expansion will be driven by the hydrodynamical equations 
without source terms. 
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Some numerical examples for asymmetric initial conditions are shown on Fig. 
I and Fig. 2. 

4-2. Asymptotic hydro without gluon wind effects 

Let us observe, that the analytic results indicate a fast transverse flow at the end 
of the inflatory period. 

An asymptotic solution of the hydrodynamical problem is given by the large 
time, large volume approximation, which implies 

(18) 
(19) 
(20) 
(21) 
(22) 

Substituting this solution into the expressions for the observables, we obtain a very 
interesting result, namely that the observables freeze-out in the sense that their 
value becomes independent from the time of the observation in the asymptotically 
late period of the expansion, when the inflation effects are already over. 

(23) 
(24) 
(25) 
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i?sidc ~ ^out ~ R\ong — * tayTa/m, (26) 

independently of the time, for asymptotically large values of the time, which implies 
that the region of homogeneity becomes spherically symmetric. Due to this reason, 
it is possible to show, that the cross-terms all vanish 

R 2 s ,o ~ R 2 ,i ~ R 2 s ,i (27) 

in the late stage of the expansion, not only in the system of ellipsoidal expansion [ 
||, but also in the frame of the observation. 

Hence the observable part of the expanding fireball becomes spherically sym- 
metric, after the inflation to large source sizes all the initial ellipsoidal geometrical 
symmetry is washed out from the local regions of homogeneity. Only the thermal 
scales [ [l5[ [l6[ ^| dominate the observable HBT radii, and these scales become 
- due to the very nature of the thermal smearing - direction independent. 

4-3. Relativistic corrections 

We find that the flow field for asymptotically large times approaches a spherically 
symmetric, Bjorken type flow field, — x^/t — 7(1, v) with v = r/t. This is 
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Fig. 1. The effect of the gluon wind on the hydro solution. We plot the evolution 
of [X, Y, Z) , the principal axis of the expanding ellipsoid, and the corresponding 
measurable HBT radius parameters, {R x , R y , Rz) for expanding fireballs in hy- 
drodynamics, with inflation and source terms in the energy equation, due to a 
gluon wind. The characteristic geometrical sizes increase exponentially with time, 
while the measurable HBT radii decrease exponentially with time. The character- 
istic time scales to reach the equality for all the HBT radii become much shorter 
than on Fig. 3, which is the case without inflation. The initial conditions are 
(X ,Y ,Z ) = (5,6,2) fin, (X,Y,Z) = (0.,0.,0.3), m = 940 MeV, T lfi = 94 MeV 
and je = jo = 0.3 c/fm. 
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Fig. 2. The effect of the gluon wind on the observables for the inflating fireballs, 
the same solution as on Fig.l. During the period of the inflation, neither the slopes 
nor the HBT radii saturate. 
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60 




Fig. 3. The time evolution of (X,Y,Z), the principal axis of the expanding 
ellipsoid, and the corresponding measurable HBT radius parameters, {R x , R v , Rz) 
for expanding fireballs in hydrodynamics, without inflation and source terms. The 
initial conditions are (X ,Y ,Z ) = (6,5,2) fm, (X ,Y ,Z ) = (0,0,0), T l = 470 
MeV to = 940 MeV. The acceleration period is over after about 10 fm/c, the 
solution tends to a coasting, ellipsoidal fireball. The HBT radii approach the same, 
constant and direction independent value, although the actual source sizes keep on 
expanding and remain direction dependent. Note the qualitative similarity to the 
results of ref . [ H) . 
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Fig. 4. The time evolution of the Hubble constants (X, Y, Z) and the measurable 
slope parameters (T x , T y , T z ) for the hydro solution shown in Fig. 3. The saturation 
of the slopes of the single particle spectra happens at the same time as that of the 
Hubble constants, which happens earlier than the saturation of the measurable HBT 
radius parameters. 
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the same flow profile as the Hubble flow of our Universe. Furthermore, due to the 
spherical symmetry of the asymptotic flow profile, the local lengths of homogeneity 
will also become spherically symmetric. Due to the large value of the geometrical 
scales, the lengths of homogeneity will measure only the thermal scales [ [l5| [l(], [| Q 
, and due to the homogeneity of the temperature and the asymptotic flow profile, 
the volume of homogeneity becomes spherical, and all the thermal length scales will 

be given by R out ~ R side ~ R\ ong ^ t \f^-- Interestingly, this thermal length 
scale will become a constant of motion in the late stages of the expansion, as the 
temperature decreases in time as T{t) = Tf(tf/t) 2 , so all the radius parameters 
Ft 

tend to t 



The relativistic corrections to this result can be estimated following the lines of 
reasoning and the conditions given in ref. [ [l^, |l6| when analyzing the connection 
between the observables and a relativistic generalization of the source functions. 
Under certain conditions, spelled out in great details in ref. [ [l5| and investigated 
numerically also in ref. [ |l^| , the relativistic form of result is 

Ft 

Rout — Rside — Rlong — Tf \ , (28) 

V m t 

where Tf = y^t 2 — r 2 at the time when the accelerationless expansion sets in, m* = 
y/ m 2 + p 2 + p 2 is the transverse mass and Tf is the value of the local temperature 
in the center of the transverse plane at the propertime 77. 

4-4- Microscopic dynamics 

It is interesting to note, that the time evolution of the observables from the ellip- 
soidally symmetric hydro solution is very much similar to the time evolution given 
by T. Humanic's cascade code, where the pions, kaons, nucleons and other hadronic 
resonances take part in rescattering, starting from a very hot and dense initial 
state [ 25 . The realization that the HBT radii may reach a saturation limit, but 
somewhat later, than the single particle spectra, was recently predicted in ref. [ p6| , 
and confirmed in Figs. 3, 4. Strong transverse push was generated also by massive 
quark matter produced in the middle of Au+Au collisions from a cascade model in 



ref. [ |27| . It seems that a various gases of massive quanta that undergoes strong 
rescattering yields qualitatively similar results to the hydrodynamic solution that 
we described here. 



5. Summary and Conclusions 

We find that an inflatory period automatically leads to a spherical symmetry for 
the observable HBT radius parameters in an analytically solvable model of fireball 
hydrodynamics with source terms motivated by the flux of high-pt gluon wind 
through the expanding hadronic matter, created in high energy collisions of heavy 
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ions. The effect is very similar in nature how a flat and spatially homogeneous 
Universe is obtained after a period of inflation in astrophysics. Our results underline 
the similarity between the physics of 'Little Bangs' at the smallest experimentally 
accessible scales and the physics of the 'Big Bang' of our Universe, at the largest 
observable scales. 

We also find, that the form of the observables is exactly the same for cases 
with or without inflation, and for various changes of the equation of state in the 
intermediate steps of the time evolution. This implies that various different equa- 
tions of state from various different initial conditions may lead to exactly the same 
hadronic final state. Furthermore, the time evolution of the parameters of the hy- 
drodynamical solution is sensitive not only to the equations of state but also to the 
magnitude of the source terms in the energy and momentum balance equations, 
which may provide non-equilibrium mechanisms to connect the initial state to a 
violently exploding final state. 

These results imply that i) the hadronic final state of fireball evolution does 
not remember the path (hence any earlier transient phase of matter) of its time 
evolution, but ii) given the equation of state, and given the possible source terms 
in the energy and the Euler equation, the initial state can still be reconstructed, 
from the observables in the final hadronic state. 

Further studies are in progress in for the publication of the ellipsoidally symmet- 
ric solutions of relativistic hydrodynamics and for the calculation of the observables 
corresponding to these relativistic generalizations of the present results. 
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